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ABSTRACT 

We construct a new family of analytic models of black hole accretion disks in dynamical equilibria. Our construction is based on 
assuming distributions of angular momentum and entropy. For a particular choice of the distribution of angular momentum, we 
calculate the shapes of equipressure surfaces. The equipressure surfaces we find are similar to those in thick, slim and thin disks, and 
to those in ADAFs. 
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1. Introduction 

In accretion disk theory one is often interested in phenomena 
that occur on a "dynamical" timescale To much shorter than the 
"viscous" timescale T[£] needed for angular momentum redis- 
tribution and the "thermal" timescale T[S] needed for entropy 
redistributiorfl, 



(1) 



The question whether it is physically legitimate to approximately 
describe the black hole accretion flows (at least in some "aver- 
aged" sense) in terms of stationary (independent on t) and axi- 
ally symmetric (independent on cp) dynamical equilibria, is not 
yet resolved. While observations seem to suggest that many real 
astrophysical sources experience periods in which this assump- 
tion is quite reasonable, several authors point out that the re- 
sults of recent numerical simulations seem to indicate that the 
MRI and other instabilities make the black hole accretion flows 
genuinely non- steady and non- symmetric, and that the very con- 
cept of the separate timescales (Q]) may be questionable in the 
sense that locally To ~ T~[X] ~ T[S\. However, this assump- 
tion has been made in all existing comparisons between the- 



ory and observations, be they by detailed spectral fitting (e.g. 
Sha fee & al]|2007U2008l and refe rences there), line profile fit- 
ting (e.g. lFabian & Vau ghan 2003), or studying small amplitude 
oscillations (see Ab ramowic zl 120051 for references). It seems 
that the present understanding of the black hole accretion phe- 
nomenon rests, in a major way, on studies of stationary and axi- 
ally symmetric models. 

From the point of view of mathematical self-consistency, in 
modeling of these stationary and axially symmetric dynamical 
equilibria, distributions of the conserved angular momentum and 
entropy, 



(2) 



* permanent address 

1 We use the spherical Boyer-Lindquist coordinates t, 0, r, the ge- 
ometrical units c = 1 = G and the + signature. The Kerr metric 

is described by the "geometrical" mass M and the "geometrical" spin 
parameter < a < 1, that relate to the "physical" mass and angular 
momentum by the rescaling, M - GM phys /c 2 , a - J vhys /(Mc). Partial 
derivatives are denoted by d t and covariant derivatives by V,-. 



may be consi dered as being free functions of the L a grangian 
coordinates (Ostriker & al. 1966; Ab ramowic 3 119701: [Bardeen 
Il970l) . The Lagrangian coordinates rj are defined by demand- 
ing that a narrow ring of matter (£, £ + d£), (rj, rj + drf) has the rest 
mass dMo = po(£, rj)d%dr] with po being the rest mass density. In 
the full physical description, the form of the functions in (|2]) is 
not arbitrary but given by the dissipative processes, like viscos- 
ity and radiative transfer. At present, several important aspects 
of these processes are still unknown, so there is still no practi- 
cal way to calculate physically consistent models of accretion 
flows from first principles, without involving some ad hoc as- 
sumptions, or neglecting some important processes. Neither the 
hydrodynamical simulations (that e.g. use the ad hoc a = const 
viscosity prescription), nor the present day MHD simulations 
(that e.g. neglect radiative transfer) could be considered satisfac- 
tory. Furthermore, the simplifications made in these simulations 



2 



Qian Lei, M.A. Abramowicz, RC. Fragile, J. Horak, M. Machida & O. Straub: The Polish doughnuts revisited 



are mathematically equivalent to guessing free functions (such 
as the entropy distribution). Bohdan Paczynski pointed out that 
it could often be more pragmatic to make a physically motivated 
guess of the final result, e.g. to guess the form of the angular 
momentum and entropy distributions. 

In practice, it is far easier to guess and use the coordinate 
distributions of the specific angular momentum and entropy, 



L 
S 



JXr,0), 
S(r,0), 



(3) 
(4) 



than the Lagrangian distributions ([2]). However, one does not 
known a priori the relation between the conserved £ and specific 
£ angular momenta (and entropy), or the functions, £ = £(r, 0), 
n = n(r, 6). Thus, assuming 0]) and is not equivalent to as- 
suming ([2]), and usually it should be a subject to some consis- 
tency conditions. We shall return to this point in Section [6j 

In several "astrophysical scenarios" one indeed guesses 
a part icular form of @]) and flU). For example, the cele- 
brated Sha kura & Sunyaevl (Il973l) thin disk model assumes the 
Keplerian distribution of angular momentum, 



£(r,0) = £ K (r) = 



M 1 ' 2 (r 2 - 2aM l l 2 r 1 ' 2 + a 2 ) 
r 3/2 - 2Mr 1 ' 2 + aM 1 ' 2 ' 



(5) 



and the popular cold-disk-plus-hot-corona model assumes a low 
entropy flat disk surrounded by high entropy, more spherical 
corona. These models contributed considerably to the under- 
standing of black-hole accretion physics. 

The mathematically simplest assumption for the angular mo- 
mentum and entropy distribution is, obviously, 



JXr, 0) = Xo = const, 
S(r, 6) = So = const. 



(6) 
(7) 



This was used by Paczyns ki and his Warsaw team to introduce 
the th ic k disk models {A bramowicz & al.ll 1978 : Kozlowski & alJ 

19801: 



Jaroszynski & al. 



Pacz ynski & Wiital 



1980; 



Abramowicz & al.lfl98ll: I Abramowicz] 1 1 9 8 ll lPaczvnskilll982h . 
Thick disks have characteristic toroidal shapes, resembling a 
doughnut. Probably for this reason, Martin Rees coined the 
name of Polish doughnut^ for them. 

Figure [T] shows a comparison of a state-of-art MHD simula- 
tion of black-hole accretion (time and azimuth averaged) with a 
Polish doughnut corresponding to a particular Xo- Both models 
show the same characteristic features of black hole accretion: 
(i) a funnel along the rotation axis, relevant for jet collimation 
and acceleration; (ii) a pressure maximum, possibly relevant for 
epicyclic oscillatory modes; and (iii) a cusp-like self-crossing of 
one particular equipressure surface, relevant for an inner bound- 
ary condition , and for stabilization of the Papal oizou-Pringle 
(Blaes 1987), thermal, and viscous instabilities (Abram owicz I 
1971). The cusp is located between the radii of marginally sta- 
ble and marginally bound circular orbits, 



Fmb < T" cusp < ^ms — IS CO. 



(8) 



Polish doughnuts have been useful in semi-analytic studies of 
the astrophysical appearan ce of super-Eddington acc retion (see 
e.g. ISikorall 19711: iMadaul 119881: ISzuszkiewicz & alJl!996b and 
in analytic calculations of small- amplitude oscilla tions of ac- 
cretion structures in connection with QPOs (see e.g. lBlaes & al.l 



2 However, real Polish doughnuts (called pyczki in Polish) have 
spherical shapes. They are definitely non-toroidal — see e.g. 



Z-Axis 




Fig. 1. Equipressure surfaces in a very simple and analytic Polish 
doughnut (left, with linear spacing), and a sophisticated, state- 
of-art full 3D MHD numerical simulation (right, with logarith- 
mic spacing). Although the shapes of equipressure surfaces are 
remarkably similar, in the numerical model the pressure gradient 
is seriously larger, and visibly enhanced along roughly conical 
surfaces, a pproximately 30° from the eq uatorial plane. (Figure 
taken from I Abramowicz & Fragilell2008l) 



http : //en. wikipedia. org/wiki/Paczki 



2006). In the same context, numerical studies of their oscillation 
properties for diffe rent angular momentu m distributions were 
first carried out by iRezzolla & al.l (l2003allbh . Moreover, Polish 
doughnuts are routinely used as convenien t starting initial con- 
figurations in numerical si mulations (e.g. Hawley et al. 20011: 
iDe Villiers & Hawlevll2003h . Recently, [Komissarov (2006) has 
constructed analytic models of magnetized Polish doughnuts. 
However, a closer inspection of Figure [T] reveals that the numer- 
ically constructed model of accretion has a (much) larger "verti- 
cal" pressure gradient than the analytic Polish doughnut, and that 
in the numerical model the gradient is visibly enhanced along 
roughly conical surfaces, approximately 30° from the equato- 
rial plane. This (and several other) detailed features of the ac- 
cretion structure cannot be modeled by either the Keplerian nor 
the constant angular momentum assumption alone. We suggest 
and discuss in this paper a simple but flexible ansatz, that is a 
combination of the two standard distributions, Keplerian © and 
constant ©. The new ansatz preserves the virtues of assuming 
the standard distributions where this is appropriate, but leads to 
a far richer variety of possible accretion structures, as are indeed 
seen in numerical simulations. 



2. Assumptions and definitions 

We assume that the accretion flow is stationary and axially 
symmetric. This assumption expressed in terms of the Boyer- 
Lindquist spherical coordinates states that the flow properties 
depend only on the radial and polar coordinates r, 6, and are in- 
dependent on time t and azimuth cp. We also assume that the 
dynamical timescale is much shorter than the thermal and vis- 



Qian Lei, M.A. Abramowicz, RC. Fragile, J. Horak, M. Machida & O. Straub: The Polish doughnuts revisited 



3 



cus ones ([]]). Accordingly, we ignore dissipation and assume the 
stress-energy tensor in the perfect fluid form, 



T\ = (p + e)u l u k - pd\ 



(9) 



with p and e being the pressure and energy density, respectively. 
The four velocity of matter u l is assumed to be purely circular, 



u l = (u\ u*, 0, 0). 



(10) 



The last assumption is not fulfilled close to the cusp (see Figure 
[U, where there is a transition from "almost circular" to almost 
"free-fall" radial trajectories. Nevertheless, the transition could 
be incorporated in the form of the inner boundary condition ( the 
relativistic Roche lobe overflow, see e.g. Abramo w iczll 1 9 8 5l) . 

One introduces the specific angular momentum £, the angu- 
lar velocity Q, and the redshift factor A by the well known and 
standard definitions, 



u t w 



(id 



The specific angular momentum and angular velocity are linked 
by 



(12) 



Qgt<j> + gtt £gt<i>+g<p<p 

The conserved angular momentum £ is given by, 

£= (P±^L £ , (13) 
Po 

3. The shapes of the equipressure surfaces 

In this section we briefly discuss one p articularly useful re- 
sult obtained by iJaroszynski & al.l ([1980). It states that for a 
perfect fluid matter rotating on circular trajectories around a 
black hole, the shapes and location of the equipressure surfaces 
p(r, 6) = const follow directly from the assumed angular mo- 
mentum distribution © alone. In particular, they are indepen- 
dent of the equation of state, p = p(e, S), and the assumed en- 
tropy distribution ©. 

For a perfect-fluid matter, the equation of motion V/ T l k = 
yields, 



dtp 



ld i g tt -2£d i g u > + J?d i g+ 



p + e 2 g tt - 2£ + £ 2 gW 

which may be transformed into, 

dtp £diQ 
= d; In A + 



1 -££l 



(14) 



(15) 



p + e 

From the second derivative commutator d r de - ded r of the above 
equation, 

d r p d e e - d e p d r e _ d r Q d e £ - d e £l d r £ 
(p + 6) 2 " (1-XQ) 2 ' 

one derives (see e.g. Abramowicz ] fl97l the von Zeipel 
condition: p(r, 6) = const surfaces coincide with those of 
e(r, 6) =const, if and only if the surfaces £(r, 6) = const coincide 
with those Q(r, 6) = consfl Obviously, the constant angular mo- 
mentum case satisfies the von Zeipel condition. 



3 The best known Newtonian version of the von Zeipel condition 
states that for a bary tropic fluid p = p(e), both angular velocity and 
angular momentum are constant on cylinders, Q = Q(R), £ 
with R = r sin 6 being the distance from the rotation axis. 



Jaroszvnski &alJdl980h have also discussed a general, non 
barytropic case. They wrote equation (ITU) twice, for i = r and 
i = 0, and divided the two equations side by side to get 



dr_p_ = d r g tt -2£d r g t * + £ 2 d r g 
d e p deg tt -2£d e g t *+£ 2 



= ~F (r, 6) . 



(17) 



For the Kerr metric components one knows the functions g lk = 
g lk (r, 0), and therefore the function F(r, 6) in the right hand side 
of (IT71) is known explicitly in terms of r and 6, if one knows or 
assumes the angular momentum distribution £ = £(r, 0). This 
has an important practical consequence. 

Let = 6(r) be the explicit equation for the equipressure 
surface p(r,0) =const. It is, d6/dr = -d r p/dep. If the function 
F(r, 6) in dT71) is known, then equation ([T71) takes the form of an 
ordinary differential equation for the equipressure surface, 6 = 
0{r\ 



d6 
dr 



(18) 



with the explicitly known right hand side. It may be therefore 
directly integrated to get all the possible locations for the equi- 
pressure surfaces. 

4. The angular momentum distribution 

4.1. Physical arguments: the radial distribution 



Jaroszvnski "&aD([T98Qh discussed general arguments showing 
that the slope of the specific angular momentum should be be- 
tween two extreme: the slope corresponding to £ = const and 
the slope corresponding to Q = const. These two cases, together 
with the Keplerian one £ = £k, may be considered as useful 
archetypes in discussing arguments relevant to the angular mo- 
mentum distribution. 

Indeed, far away from the black hole r » rp , these argu- 
ments are well known (see e.g. [Frank TiD[2002) and together 
with numerous numerical simulations show that typically (i.e. in 
a stationary case with no shocks) the specific angular momen- 
tum should be slightly sub-Keplerian £(r,n/2) « £k(t). There 
is a solid consensus on this point. 

The situation close to the black hole is less clear because 
there is not sufficient knowledge of the nature of the stress op- 
erating in the innermost part of the flow, i.e. approximately be- 
tween the horizon and the ISCO. Formally, one may consider 
two extreme ideal cases, depending whether the stress is very 
small or very large. 

In the first case, the almost vanishing stress implies that the 
fluid is almost free-falling, and therefore the angular momentum 
is almost constant along fluid lines. This l eads to £(r,n/2) » 
const. Such situa tion is typical for the th in IShakura & Sunyaevl 
(Il973l) and slim (I Abramo wicz & al .II 1 9 8 8l) accretion disks. In the 
second case, one may imagine a powerful instability like MRI, 
which occurs when d£l/dr ± 0. It may force the fluid closer to 
the marginally st able state Q = const. This situation may be rele - 
vant for ADAFs dNaravan & Yilll 98 8U Abramowicz & alJI 1995b . 

4.2. The new ansatz 

We suggest adopting the following assumption for the angular 
momentum distribution, 



£(<R\ JXr, 6) 



£o(^f sin 2 ^ for r>r m / 



£ m ,(r)sin 2 ^ 



for r < r m 



(19) 
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Fig. 2. (a) and (b): the distribution of angular momentum on the equatorial plane. Thick lines correspond to the angular momentum 
predicted by our analytic formula, dashed l ines show t he Ke plerian angular momentum distribution and dot s to the simulation data, 
(a): Ke rr geometry a = 0.9 simulations by lSadowskil (|2008|) . (b): Pseudo-Newtonian MHD simulations by iMachida & Matsumotol 
( 2008). (c) and (d): angular momentum off the equatorial plane, normalized to its equatorial plane value, A = X(r, 6) /X(r, n/2). Lines 
correspond to the sin 2r distribution at r = 10r G : long-dashed y = 0.5, dotted y = 1.0, and short dashed y = 1.5. Points are taken 
from time-dependent, fully 3-D, MHD numerical simul ations. They correspond to time and azimuthal averages at the same radial 
location, r = lOrc. (c): Points from the simulations of Mach ida & Matsumotol (|2008l) in the Paczynski- Wiita po tential — triangle: 
High temperature case; square: Low temperature case, (d): Points from the simulations of iFragile & al.1 (120071) in Schwarzschild 
(dots) and a = 0.9 Kerr (crosses) spacetimes. 



The constant Xo is defined by Xo = ^7-C^(r m5 ). For the "hydro- 
dynamical" case, the function £ ms (r) is constant, 

£ ms (r) = Xo UK(r ms )/£of = const, (20) 

while for the "MHD" case its is calculated from the £l(r) = 
^KiTms) = const condition, 

&ms gw( r > 71 V 2 ) + 8t<p(r, n/2) 



£ms(r) = — 



(21) 



&ms gt<f>(r, n/2) + g tt (r, tt/2) 

Thus, there are only three dimensionless parameters in the 
model: (J3, y, rf). Their ranges are, 

0<)8<1, -l<y<l, l</7<^. (22) 

The function JLk(v) is the Keplerian angular momentum in the 
equatorial plane, 6 = n/2, which for the Kerr metric is described 
by formula © and rj max = XH^)/ '£,K(r ms ). An equipressure 
surface that starts from the cusp is marginally bound for /3 = 0, 
y = and 77 = r] max . 

4.3. Angular momentum on the equatorial plane 

On the equatorial plane, sin = 1 , and therefore only ft and 77 
(through Xo) enter the distribution formulae ([T9b . 



JXr,n/2) = 



Jin 



for r <r m 



(23) 



When p = 0, the angular momentum is constant, £ = Xo, and 
when p = 1, it equals the Keplerian one, X = JLk- 

For small values of J3 the assumed equatorial plane angular 
momentum (|23l) reproduces the characteristic shape, shown in 
Figure [2 which has been found in many numerical simulations 
of accretion flows — including stationary, a xially symmetric, a 
visco sity, hydrodynamical "slim disks" (e.g. Abr amowicz & al.l 
Il988l) . an d more recent, fully 3-D, non- s t ationary MHD simu la- 
tions (e.g. IMachida & Matsumotol 120081: IFragile & al.ll2008l) . It 
corresponds to a distribution that is slightly sub-Keplerian for 
large radii, and closer to the black hole it crosses the Keplerian 
distribution twice, at r center > r ms and at r cusp < r ms , forming a 
super- Keplerian part around r ms . For r < r cusp the angular mo- 
mentum is almost constant. 



4.4. Angular momentum off the equatorial plane 

Numerical simulations show that away from the equatorial 
plane, the angular momentum fall s off. Figure 12 shows that in- 
deed several MH D simulations (iMachida & Matsumoto 2008; 
Fra gile &~aT1l2008L Figure 2c and 2d respectively), feature a drop 
of angular momentum away from the equatorial plane. This be- 
havior is reflected by the term sin 2y # in (|T9l ). One may see 
that this form a ccurately mimics the outco me of the numeri- 
cal simulations. iProga & Begelmanl (l2003allbb also studied ax- 
isymmetric accretion flows with low specific angular momen- 
tum usin g numerical simulat i ons. In their inviscid hydrodynam- 
ical case Proga & Begelman (2003a) found that the inner accre- 
tion flow settles into a pressure-rotation supported torus in the 
equatorial region and a nearly radial inflow in the polar funnels. 
Furthermore, the specific angular momentum in the equatorial 
torus was nearly constant. This beh avior changes once magneti c 
fields are introduced, as shown in IProga & Begelma n ( 2003b). 
In the MHD case, the magnetic fields transport specific angular 
momentum so that in the innermost part of the flow, rotation is 
sub-Keplerian, whereas in the outer part, it is nearly Keplerian. 
Similar rotational profiles are also found in MHD simulations 
of the collapsar m odel of gamma-ray bursts (Pro ga & a fl l2003l: 
iBaiotti &~al]|2008l) . which use a sophisticated equation of state 
and neutrino cooling (instead of a simple adiabatic equation of 
state). Therefore, it appears that the rotational profile assumed in 
our model is quite robust as it has been obtained in a number of 
numerical experiments with various microphysics. 



5. Results 

Figures [3] and H] show sequences of models calculated with the 
new ansatz (IT9l) for black-hole spins a - and 0.5, respectively. 
For these models we hold rj = r] max fixed, while j3 and y are 
varied over the limits of their accessible ranges. 

5.1. Equipressure surfaces on the axis of rotation 

Figures [3] and 01 show an interesting change of the behavior of 
equipressure surfaces close to the axis with increasing y. No 
equipressure surface can cross the symmetry axis when the de- 
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Fig. 3. Equipressure surfaces for a = and rj = rj max = 1.085. Five rows correspond to ft = (0.0), (0.1), (0.5), (0.9), (0.99) from the 
top to the bottom. Four columns correspond to y = (0.0), (0.1), (0.5), (0.9) from the left to the right. The upper left corner shows a 
"standard" Polish doughnut. The lower right corner shows an almost Keplerian disk at the equatorial plane, surrendered by a very 
low angular momentum envelope. 



pendence of the angular momentum on 6 is weak. This is the 
case for the first three columns of Figure[3]where y < 0.5. On the 
other hand, for the angular momentum distributions with higher 
y the equipressure surfaces cross the axis perpendicularly. This 
happens in plots of the last column of Figure [3] This behavior 
can be understood easily from the limit of rd6/dr as 6 — > 0. In 
Schwarzschild spacetime equations (TT71) and (fTSt give 



rdO 
lim 

e^o + dr 



2£ 2 K (r) 



lim (sm 4y - 3 0). 



(24) 



The limit on the right-hand side is either 0, 1 or oo, depending on 
the value of y. When y < 3/4, rd6/dr = and no equipressure 
surface goes across the axis. On the other hand, when y > 3/4 
equipressure surfaces cross the axis perpendicularly. Of course, 
a stationary torus may exist only within an equipotential sur- 
face located inside the Roche lobe, i.e., the critica l self-crossing 
equipotential within the cusp (I Abramowicz II 1 9 8 5b . 

5.2. Comparison with numerical simulations 

Figure[5]illustrates that the results of the analytic models are well 
matched with results of modern 3-D MHD numerical simula- 



tions (here taken from Fragile & al. 2007, 2008). For the correct 
choice of parameters, the model can reproduce many of the rele- 
vant features of the numerical results, including the locations of 
the cusp and pressure maximum, as well as the vertical thickness 
of the disk. At this stage, such qualitative agreement is all that 
can be hoped for. One notable difference between the analytic 
and numerical solutions is the behavior inside the cusp. While 
the analytic equipressure surfaces formally diverge toward the 
poles, the numerical solution maintains a fairly constant verti- 
cal height, which is also evident in Figure [T] This is because in 
the region inside the cusp, our assumption (ITOl) about the form 
of the velocity is not valid — velocity cannot be consistent with 
a pure rotation only, u l = (u\ u? 9 0, 0). In this region the radial 
velocity u r must be non-zero and large. Thus, accuracy of our 
analytic models may only be trusted in the region outside the 
cusp, r > r cusp . 



6. Discussion 

In this paper we assumed a form of the angular momentum 
distribution (|T9l) and from this calculated the shapes and loca- 
tions of the equipressure surfaces. This may be used in calcu- 




Fig. 4. Equipressure surfaces for a = 0.5 and rj = rj max = 1.079. Five rows correspond to J3 = (0.0), (0.1), (0.5), (0.9), (0.99) from the 
top to the bottom. Four columns correspond to y = (0.0), (0.1), (0.5), (0.9) from the left to the right. The upper left corner shows a 
"standard" Polish doughnut. The lower right corner shows an almost Keplerian disk at the equatorial plane, surrendered by a very 
low angular momentum envelope. 



lating spectra (in the optical ly thick case) by t he same "surf ace" 
method as used in works bv lSikoral (Il97ll) and lMadaul (119881) . 

We plan to construct the complete physical model of the in- 
terior in the second paper of this series. Here, we only outline 
the method by considering a simplified toy model. Let us denote 
p = 6 + p. We assume a toy (non-bary tropic) equation of state 
and an entropy distribution, by writing, 



K = K(r, 6). 



(25) 



p = e K ^p, 

Let us, in addition, define two functions connected to the entropy 
distribution, 



d e K = K(r,0). ^-=A(r,6), 
d e K 



(26) 



From the obvious condition that the second derivative commuta- 
tor of pressure vanishes, (d r d e - d e d r )p = 0, and equations d25lh 
(f26l) and (IT71) one derives, 

d r Gq - do G r 



k = -- 

G r — A Go 

where G r and Go are defined as 

G,-(r, 6)= — 
P 



(27) 



(28) 



and can be calculated from the angular momentum distribution 
using equation (fT4l) . From (f27l) it is obvious that one cannot in- 
dependently assume the functions K(r, 6) and A(r, Assuming 
A(r, 6) is equivalent with assuming the shapes of isentropic sur- 
faces. Indeed, from (f26t one concludes that the function 6 = 
6s(r) that describes an isentropic surface is given by the equa- 
tion, 



dO 
dr 



-A(r,6). 



(29) 



that may be directly integrated. Then the condition (l27l) gives 
the physical spacing ("labels") to the isentropic surfaces, and 
through the equation of state (f25l) also to equipressure surfaces 
and isopicnic (p = const) surfaces. 

Note, that a possible choice A = G r /Go corresponds, obvi- 
ously, to the "von Zeipel" case in which equipressure and isen- 
tropic surfaces coincide. In this case the denominator in (f27l) 



4 A somewhat similar si tuation in the case of rotating stars is known 
as the von Zeipel paradox (Tassoul 1978): Pseudo-barytropic models in 
a state of permanent rotation cannot be used to describe rotating stars 
in strict radiative equilibrium. 
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Fig. 5. Comparison of pressure distributions between the an- 
alytic model (dark lines) and numerical si mulations (colors) . 
The r esults of MHD simulations (taken from Frag ile & al.ll2007i 
2008) have been time-averaged over one orbital period at r = 
25rc. Upper panel: Schwarzschild black hole (a = 0); the ana- 
lytic model parameters are n = 1.085, ft = 0.9, and y = 0.18. 
Lower panel: Kerr black hole (a = 0.5); the analytic model pa- 
rameters are n = 1.079, j3 = 0.7, and y = 0.2. 

vanishes, implying a singularity unless the numerator also van- 
ishes. The condition for the numerator to vanish is, however, 
equivalent to the von Zeipel condition. 

7. Conclusions 

The new ansatz dT9b captures two essential features of the angu- 
lar momentum distribution in black hole accretion disks: 

1. On the equatorial plane and far from the black hole, the 
angular momentum in the disk differs only little from the 
Keplerian one being slightly sub-Keplerian, but closer in it 
becomes (slightly) super- Keplerian and still closer, in the 
plunging region, sub-Keplerian again and nearly constant. 

2. Angular momentum may significantly decrease off the equa- 
torial plane, and become very low (even close to zero, in a 
non-rotating "corona"). 

Models of tori described here may be useful not only for accre- 
tion disks but also for tori that form in the latest stages of neu- 
tron star bi nary mergers. This is relevant f or gamma ray bur sts 
(IWitt&alJI 19941) and gravitational waves dBaiotti & al.ll2008h . 
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